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We report a detailed first-principles study of the structural, electronic and vibrational properties of 
the superconducting C32 phase of the ternary alloy CaAb-zSiz, both in the experimental range 0.6 < 
x < 1.2, for which the alloy has been synthesised, and in the theoretical limits of high aluminium 
and high silicon concentration. Our results indicate that, in the experimental range, the dependence 
of the electronic bands on composition is well described by a rigid-band model, which breaks down 
outside this range. Such a breakdown, in the (theoretical) limit of high aluminium concentration, is 
connected to the appearance of vibrational instabilities, and results in important differences between 
CaAl2 and MgB2. Unlike MgB2, the interlayer band and the out-of-plane phonons play a major 
role on the stability and superconductivity of CaAISi and related C32 intermetallic compounds. 

PACS numbers: 74.70.Dd, 71.15.Nc. 74.25.Kc, 74.25.Jb 
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The discovery^ of superconductivity at T c = 39 K in 
magnesium diboride (MgB 2 ) is remarkable not only for 
its high critical temperature relative to conventional su- 
perconductors, but also because it has destroyed preju- 
dices and created new hopes on electron-phonon (e-ph) 
superconductors. MgB 2 crystallises in the AlB 2 -type 
structure, also called C32 phase (P6/mmm, No. 191), 2 
forming hexagonal layers of graphene-like sheets of B 
atoms intercalated with Mg planes; the synthesis and 
study of materials with the A1B 2 structure has therefore 
attracted a great deal of interest in the last few years. 

In the present paper we will focus on the ternary al- 
loy CaAl^zSia;, which has been recently synthesised 3 - for 
a relatively wide range of composition 0.6 < x < 1.2 
and, in this range, has been reported to be isostruc- 
tural to MgB 2 , with honeycomb, graphene-like sheets of 
randomly distributed silicon and aluminium atoms, and 
alkaline-earth atoms intercalated between these sheets^ 
Superconductivity has been detected in the whole com- 
position range; the critical temperature T c (x) strongly 
depends on composition, showing a sharp peak at x = 1 
(where T c = 7.8 K), and dropping off rapidly as the Al/Si 
ratio departs from unity. 3 

Several experiments^*^ and ab initio calcula- 
tions£*£**^fii have concentrated on the x = 1 case (i.e. 
CaAISi). Concerning the superconducting mechanism, 
Mazin et al& have suggested that the available exper- 
imental data could be reconciled with a traditional 
e-ph scenario, assuming the existence of a soft phonon 
mode (~ 30-40 cm -1 ); recently Huang et alm^ have 
calculated the phonon spectrum and electron-phonon 
coupling of CaAISi, revealing the existence of a very soft 
(and, at some wavevector, even unstable) B\ g branch, 
which plays an important role in its superconducting 
properties; they also claim that this mode, which gets 
softer under pressure, can be invoked to explain the 



In contrast to the x = 1 case, neither the wide stability 
range 0.6 < x < 1.2 of of the C3 2 phase, nor its insta- 
bility outside this range (the limit of high aluminium or 
silicon concentration) , nor any other physical property of 
CaAl 2 _ x Si x has been, so far, theoretically studied as a 
function of composition, at least to our knowledge. 

Our paper reports the first detailed ab initio analysis 
of the structural, electronic and vibrational properties 
of this ternary alloy, both in the composition range for 
which it has been actually synthesised, and outside it, in 
the theoretical limits of high aluminium and high silicon 
concentration. The former limit (x — > 0) is of particular 
interest, since Ca and Al have the same valence config- 
uration as Mg and B, respectively, and, therefore, the 
hypothetical*^! C32 phase of CaAl 2 would not only be 
isostructural, but also isoelectronic to MgB 2 , attracting 
attention as a possible new superconductor or, at least, 
as a relevant contribution to the understanding of MgB 2 
and related superconducting alloys^ 

The work is organised as follows. In Sec. [I] we briefly 
describe the computational details of our study, includ- 
ing the virtual crystal approximation adopted for a dis- 
ordered alloy. In Sec. [H] we discuss the electronic bands 
of the x=l case, considering both an ordered phase with 
in-plane Al-Si ordering and the disordered alloy, and com- 
pare our results with previous calculations*^ and ARPES 
data£ We also introduce a third, fictitious compound 
□ 2+ AlSi, with calcium ions replaced by a uniform back- 
ground of positive charge, whose electronic bands will 
help us identify the role of calcium ions. The vibrational 
properties and the electron-phonon coupling in CaAISi 
and D 2+ AlSi are presented in Sec. IIIII while, in Sec. IIVI 
we describe how the different Al/Si content modifies the 
electronic and structural properties of CaA^-^SLc. 
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I. COMPUTATIONAL DETAILS 

The electronic, structural and vibrational properties of 
all our crystals were calculated in the framework of the 
density functional theoryiStil and of the density func- 
tional perturbation theory^, with Kohn-Sham wavefunc- 
tions expanded on a plane- wave basis limited by a 40 Ryd 
cutoff energy42i211 We used Troullier-MartinaSl norm- 
conserving pseudopotentials generated by the FHI98PP 
package^ 2 , and consistently employed the GGA-PBE 
functional^ 3 - to approximate the exchange-correlation en- 
ergy functional. The k-space integration (electrons) was 
approximated with al2xl2xl2 Monkhorst-Pack grid 24 
for the self-consistent cycles; with the more accurate 
tetrahedron method 2 ^ and a 20 x 20 x 20 mesh for the elec- 
tronic density of states (DOS); a 18 x 18 x 18 Monkhorst- 
Pack grid for the electron-phonon coupling. After setting 
the smearing parameter— to 0.03 Ryd, the electronic to- 
tal energies were converged to 0.1 mRy, while the phonon 
frequencies were converged within 5 cm -1 with respect to 
the k-point sampling; we checked that a different choice 
of the smearing parameter between 0.02 Ryd and 0.04 
Ryd would not change our frequencies by more than 5 
cm -1 . Dynamical matrices were evaluated for phonon 
wave vectors qona6x6x6 grid, from which the phonon 
dispersion was obtained by Fourier interpolation. 

To describe a disordered CaAla— x Sia; alloy we em- 
ployed the Virtual Crystal Approximation (VCA) 2 ^, thus 
considering a P6 / mmm structure where the sites of the 
hexagonal graphene sheets are occupied by an alchem- 
ical element whose ionic pseudopotential is defined as 
v£n = (1 - x/2)v^ + (x/2)vfj n . The use of the VCA 
permits us to investigate the whole composition range 
x £ [0,2] without introducing superstructures (which on 
the other hand were not experimentally revealed in the 
samples examined with a powder x— ray diffraction tech- 
nique)^ The VCA should be justified in our case, since 
Al and Si are neighbours in the periodic table and have 
a very similar electronic structure. 

II. ELECTRONS IN CaAISi, Ca(AlSi) AND D 2+ AlSi 

The x = 1 case requires a thorough discussion, since 
this is the composition to which most experimental and 
theoretical results refer. Furthermore, there is still an 
open question concerning the actual crystal structure, 
whether a chemically disordered phase ( where each site 
of the hexagonal layers is randomly occupied by either 
Al or Si, with equal probability), or a crystal structure 
with some Al/Si ordering. The sharp x-ray diffraction 
patterns and the relatively high phase purity measured 
in CaAISi (as opposed to the case of CaA^-^Sia; with 
x^= 1) suggest a possible Al/Si ordering for x — l. 3 This 
is also what Mazin et al& suggest, proposing a crystalline 
structure with atomic in-plane ordering (i.e., graphene- 
like sheets with alternated Al and Si atoms) of P6to2 
(No. 187) symmetry^ obviously lower than the C32 sym- 



metry described above, and apparently favourable with 
respect to other possible ordering patterns. 

An analysis of powder x-ray diffraction patterns per- 
formed by Imai et al£§> suggests, instead, the disordered 
phase. The question is still under debate and no conclu- 
sive experimental measure has been reported yet. 

To elucidate this controversy we have repeated Mazin's 
first-principles calculations, considering both the ordered 
P6m2 structure and the disordered P6/mmm VCA 
phase. Our comparison of these results with the ARPES 
experiments^ gives additional arguments in favour of the 
ordered phase. Moreover, the VCA phase represents a 
crucial reference for our subsequent rigid-band model. 
Hereafter we will indicate our disordered VCA phase as 
Ca(AlSi) and the ordered phase as CaAISi. The CaAISi 
and Ca(AlSi) lattice parameters were obtained by opti- 
mizing the structural parameters c and a up to a maxi- 
mum residual stress of 1 MPa; the numerical values are 
listed in Table [Q 

TABLE I: Experimental lattice parameters at NTP^ and op- 
timised lattice parameters for the ordered CaAISi and the 
VCA-disordered Ca(AlSi) obtained in this work, based on 
GGA-PBE exchange-correlation and pseudopotentials. 





a(a.u.) 


c (a.u.) 


c/a 


Exp 


7.916 


8.315 


1.050 


CaAISi 


7.979 


8.446 


1.058 


% error 


0.79 


1.57 


0.76 


Ca(AlSi) 


7.965 


8.485 


1.065 


% error 


0.62 


2.04 


1.42 



The agreement of the calculated lattice parameters 
with the experiment is good (slightly better for the or- 
dered phase), with errors of the order of 1%. The cor- 
responding electronic bands are shown in Fig. ^ for the 
ordered structure (upper panel) and the disordered VCA 
phase (middle panel); they agree quite well with the pre- 
vious calculations^^, and, in particular, with the corre- 
sponding full-potential results—. 

We will start our discussion with the energy bands of 
the ordered CaAISi structure. 



CaAISi 

In CaAISi the sp 2 -hybrids of aluminium and silicon 
form three bonding a bands: the lowest occupied a band 
is formed by s-like states and is separated by a gap from 
the other two a(2p x ^ y ) bands, which are mainly formed by 
the 2p x , 2p v states of Al and Si. The a(2p x , y ) bands have 



3 



a quasi two-dimensional (2D) character, with a small dis- 
persion along k z (the T-A direction) very similar to what 
is observed in MgB2i22i2£i We will refer to the a(2p x<y ) 
bands as a bands in what follows. The 2p z states of 
Al and Si form a bonding (tt) and an anti-bonding (tt*) 
band, which are completely three-dimensional (3D). 

Because of the reduced crystal symmetry, the tt and 
tt* bands do not cross at K, but are separated by a gap 
in the spectrum. And, unlike MgB2, both the bonding tt 
and the bonding a bands are fully occupied, because here 
there is one more electron per unit cell; so these bands 
cannot play any role for superconductivity. Moreover, 
because of the low c/a ratio, the tt and tt* bands acquire 
a substantial dispersion along k z , much higher than in 
graphite and MgB2; the final result is that, in the "upper 
half of the Brillouin zone (k z ~±tt/c, A-L-H-A 
path), a good fraction of the tt* band is occupied too, 
and the antibondig tt* band, of 3D character, crosses the 
Fermi level. 

Then, as we see in the upper panel of Fig. ^ a sec- 
ond band crosses the Fermi level. It's a 3D band whose 
wavefunctions extend in the interlayer region. Its bottom 
has a prevailing Ca3<i z 2_ r 2 character, yet the band as a 
whole may be easily traced back to the so-called graphite 
interlayer band]Q therefore, in what follows, we will stick 
to the name. This band is empty in graphite and MgB2, 
but here, in CaAISi, it's partially full, because of the dif- 
ferent electron count. 

As we see in the upper panel of Fig. |3 the CaAISi 
Fermi surface resulting from these two bands has two 
topologically disconnected sheets. The first, outer sheet 
is a sixfold gear-like surface: a wide empty cylinder, 
whose axis is along k z , connected to the neighbouring 
Brillouin zones in the k xy plane along the T-M and equiv- 
alent directions by six sleeves, stretching out from its side 
surface. This sheet of the Fermi surface actually derives 
from both bands: near k z =0 its character is mainly in- 
terlayer; near k 2 =±7r/c it is mainly tt* . The second, 
inner sheet is a hexagonal lens-like structure centered at 
r. It exclusively arises from the interlayer band, which 
cuts cf (upper panel of Fig. ^) around the centre of the 
Brillouin zone (BZ). This structure is also observed in 
the ARPES data, 6 which have revealed the presence of 
occupied states around T forming a hexagonal structure. 



Ca(AlSi) 

Comparing the Ca(AlSi) and CaAISi bands of Fig. ^ 
we can notice that, because of the increased symmetry, 
the tt — tt* energy gap disappears in the VCA: now the 
tt bands cross at K, forming an electronic pocket at Ei?. 

The most important difference is that the inner sheet 
of the Fermi surface, a large hexagonal lens in the or- 
dered CaAISi, shrinks to almost nothing in the VCA (a 
tiny circular lens, see lower panel of Fig.[2J), because now 
the Fermi level cuts the interlayer band almost at its bot- 
tom, at r (see middle panel of Fig. QJ. Otherwise, the 



outer sheet is almost identical, differing from the corre- 
sponding one in CaAISi only near the K point because of 
the additional small electronic pocket. 

The lack of a sizable hexagonal lens-like structure 
makes the Fermi surface of the disordered structure less 
compatible with the experimental ARPES map. Fur- 
thermore, unlike the ordered phase, the band structure 
of the disordered Ca(AlSi) (middle panel of Fig. ^) can- 
not explain^ an ARPES high-intensity region located at 
~ —0.8 eV, derived from tt states, which is clearly visible 
in Fig. 2 of Ref. 0. Our results seem thus to confirm the 
existence of in-plane ordering in the compound synthc- 
sised with equal amounts of Al and Si. 

While the bands of the ordered phase appear in better 
agreement with the ARPES experiment at x—1, the VCA 
band structure of the disordered Ca(AlSi) , middle panel 
of Fig. represents a very convenient reference for a 
rigid-band description of CaA^-zSiz at x ^ 1, a model 
which will prove quite accurate in the experimental range 
of compositions x G [0.6; 1.2]. In a rigid-band framework 
we find it useful to report, on the right side of Fig.^ the 
Integrated Density of States (IDOS), with black dashed 
lines corresponding to the band edges of the interlayer, 
a, and tt bands. 

In this way we may see that, below 8.8 electrons, the tt 
bonding band in the k z = begins to be emptied; below 
8.5 electrons, the interlayer band is completely empty, 
and below 8.25 electrons the a bonding band in the k z = 
±7r/c plane begins to be emptied (see Section llV|) . 

a 2+ AiSi 

If a few differences in the electronic bands have helped 
us discriminate between the ordered and disordered 
phase, their similarity is overwhelming; the band struc- 
tures in the upper and middle panels of Fig.Q]are, in fact, 
almost identical. Now we want to point out that the same 
striking similarity extends to the lower panel, where we 
show the bands of n 2+ AlSi. This fictitious compound 
has the same number of electrons and the same lattice 
parameters as CaAISi, but the Ca 2+ ion has been re- 
placed by a homogeneous positive background \3 2+ . 

We see that the absence of the calcium ion doesn't 
dramatically change the band dispersion at : the tt* is 
unchanged, while the interlayer band shows slight modifi- 
cations. This indirectly confirms the interlayer character 
of this band. The important differences with respect to 
CaAISi are located above cf, where the bunch of empty 
Ca-d bands (and the corresponding large hump in the 
DOS) is clearly (and obviously) absent. 



III. PHONONS AND E-PH COUPLING IN 

CaAISi (x = 1) 

Having become convinced that an atomic in-plane or- 
dering takes place for x = 1, we will discuss the phonon 
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FIG. 1: Upper panel: energy bands (black), DOS (blue) and 
integrated DOS (IDOS, red) for CaAlSi. Middle panel: same 
as above for the disordered Ca(AlSi); the three black dashed 
lines pinpoint, along the IDOS, the number of electrons at 
which, in a rigid band scheme, a few relevant bands begin 
to be emptied or filled (see text). Lower panel: same as the 
upper panel, for the fictitious D 2+ AlSi. 



dispersion only in the ordered phase, shown in the upper 
panel of Fig- EI At the T point, in addition to the acousti- 




FIG. 2: Fermi surface and Brillouin zone of CaAlSi in the 
in-plane ordered (upper panels) and disordered (lower pan- 
els) structure. The outer sheet originates from both the 7r* 
(around A) and the interlayer band (around F). The inner 
sheet (a lens-like structure) exclusively originates from the 
interlayer band at V, and is much larger, with more evident 
hexagonal symmetry, in the ordered structure (upper panel). 
We have chosen the blue color for the outer surface of the 
gear-like structure, and the red color for its inner surface. 
The surface of the inner, lens-like structure, is also red. With 
this choice the top view (right panels) happens to show the 
7r* part of the Fermi surface in blue and the interlayer part 
in red. 



cal modes, there are four distinct optical modes: one dou- 
bly degenerate Ei u mode (vibrations of Ca atoms against 
the Al-Si planes along x, y), one A2« mode (Ca vibrations 
against Al-Si planes along z), one doubly degenerate E2 3 
mode (in-plane Al-Si bond-stretching mode), and one Bi s 
mode (Al and Si displaced along z in opposite directions). 
Table HTl compares the frequencies uj of the optical modes 
at r obtained in this work with those reported in pre- 
vious full-potential calculations. For comparison, in the 
same table we also report the corresponding phonon fre- 
quencies for MgB^&. 

The most striking feature is the presence of a very soft 
B\g mode. First intuition, inspired by the correspond- 
ing mode in MgB 2 and other diboride o??^ 4 , would place 
this mode, which only involves vertical displacements of 
the lighter Al and Si atoms, at frequencies higher than 
those involving Ca too. We find it, instead, as the lowest- 
frequency optical mode not only at the Y point (Table[nJ , 
but also in the entire "upper half of the BZ (Fig. |3J. 
This finding can be reconciled with the mass argument 
only by a much weaker vertical force constant; and in- 
deed, unlike the other compounds, CaAlSi has a partial 
filling of the ir* bands in the "upper half" of the Brillouin 
zone (k z ~ ±tt/c, A — L — H — A path), which weakens 
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TABLE II: Calculated frequencies (cm - ) of the zone-center 
optical phonons in ordered CaAISi, as obtained here (asterisk) 
and in previous calculations. MgB2 is shown for comparison. 

Bi 9 (r) Ei„(r) A 2u (r) E 2g (r) 



CaAISi* 


105 


183 


221 


444 


CaAlS* 


99 


194 


222 


455 


CaAISi 8 


100 


187 


212 


456 



MgBa** 692 335 401 585 



the bond in the vertical direction. 

Our phonon frequencies at T (Table |Hj are in good 
agreement with all previous calculations; the same agree- 
ment is found with the entire phonon spectrum ob- 
tained by Huang et al£, with a notable exception: our 
B\g branch (upper panel of Fig. |3J) has a frequency of 
~ 50 cm -1 and is almost dispersionless in the q z = ±ir/c 
plane; Huang et al. find, instead, imaginary frequencies 
along the A-L path. We find no such dynamical instabil- 
ity either there, or, in fact, in any other phonon branch 
and wavevector. Our results agree with the overwhelm- 
ing experimental evidence^&iS in favour of the structural 
stability of CaAISi, and match the remarkable guess of 
Mazin et al, that, outside T, a soft mode must be there 
at about 35-40 cm -1 for a consistent explanation of both 
normal and superconducting properties. 

In the upper panel of Fig. right side, the calculated 
Eliashberg function a 2 F(uj), shown in red, is compared to 
the phonon DOS, shown in blue. The Eliashberg function 
evidently does not follow the phonon DOS, and shows a 
pronounced peak at ~ 50 cm -1 , clearly arising from the 
(anomalously low) B± g branch. Electrons mainly couple 
with this branch; the largest contribution to the inte- 
grated superconducting parameter A will appear along 
the T-A-L path. This coupling is mainly interband in 
character: a B lg displacement mixes the interlayer and 
7T* bands, while, by symmetry, the corresponding 7r-7r 
coupling remains negligible 35,36 . 

Since, as discussed in Sec. [H] Fermi-level interlayer 
electrons are near k z — 0, while Fermi-level 7r* elec- 
trons are near k z — ±7r/c, the electron-phonon coupling 
happens at wavevectors in the "upper half of the BZ 
(q z ~ ±7r/c, A — L — H — A path). The total A, obtained 
as a weighted sum over individual phonon linewidths, is 
0.66. The critical temperature can then be estimated 
from the Dynes equation 

rp ^log ( 1.04(1 + A) \ 

Tc = T^ eXp ( v ' A- / ^(l-0.62A) J ; (1) 



In the lower panel of Fig.|3|we report the same vibra- 
tional properties (phonon dispersion, DOS and Eliash- 
berg function) as in the upper panel, but calculated for 
the fictitious system D 2+ AlSi. Here we obviously have 
only 6 out of 9 phonon branches, since the Ca atoms are 
missing. Let us focus on the three optical modes, which 
correspond to the two E2 ff branches (degenerate at T) 
and the Bi s branch of the original compound. Our re- 
sults show that, while the two Fi2 g modes are more or less 
unchanged, the substitution of Ca ions with a uniform 
background of positive charge drives the Bi s branch to- 
wards imaginary frequencies (shown as negative in Fig.[3J) 
across the entire Brillouin zone, indicating that our ficti- 
tious compound is not stable in this structure. 

The Eliashberg function a 2 F(u) (red line) is thus well 
defined for all branches except this one. At imaginary 
frequencies we may extend, as done here, its definition 
by taking their modulus. Such an extended Eliashberg 
function cannot be integrated to yield a total A, but may 
be used to illustrate the frequency distribution of the e- 
ph interaction^ which turns out to be very similar to 
the real CaAISi (upper panel). This is not surprising, 
because in D 2+ AlSi the relevant, gear-like structure of 
the Fermi surface (not shown) is very similar to that of 
ordered and disordered CaAISi (see Fig. 0. The only 
difference is that, in the fictitious D 2+ AlSi (lower panel), 
the e-ph peak corresponding to the B lg branch, is shifted 
to negative frequencies; we may therefore conclude that, 
although the Ca ion is not so important to define the 
CaAISi electronic bands at ep, it plays a crucial role for 
its structural stability, because, with respect to the equiv- 
alent jellium-like charge, a positive ion placed between 
Al-Si layers is more effective both in binding them and 
in hindering the Al and Si vertical displacements . 

In summary, we find an entire phonon branch of imagi- 
nary frequencies in the fictitious compound □ 2+ AlSi, but 
we do not find any imaginary frequency in the phonon 
spectrum of the actual (ordered) CaAISi, in contrast to 
previous, full-potential calculations. In CaAISi we only 
find a very soft branch, which suggests the proximity of a 
structural transition. This single phonon mode, which is 
soft in the "upper half of the Brillouin zone, dominates 
the electron-phonon coupling. 

We suspect that a similar situation also holds at com- 
positions other than x=l, and, in the following section, 
we explore in some detail the relation between lattice sta- 
bility and band structure in the whole composition range. 
As we will see, the results confirm our hint, at least in 
the limit of high-aluminium concentration, and provide 
a useful reference and a strong motivation for a future, 
complete study of the electron-phonon interaction as a 
function of x. 



IV. PHYSICAL PROPERTIES OF CaAla-xSi 



with a calculated logarithmic average frequency of 151 K, 
and assuming /i* = 0.1, we obtain a T c of 5.8 K, in 
reasonable agreement with the experimental T c of 7.9 K. 



In this section we will give a quantitative description 
and attempt a qualitative interpretation of how the com- 
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FIG. 3: Upper panel: Phonon dispersion, phonon DOS (blue) 
and Eliashberg function (red) a 2 F(uj) for CaAlSi. Notice that 
the Eliashberg function has a large peak around 50 cm" , 
corresponding to the soft Bi 9 branch. Lower panel: same as 
above but in the fictitious □ 2+ AlSi; the Bi 9 branch has an 
equally large e-ph coupling, but now its frequency are imag- 
inary (negative in the plot) in the entire BZ, revealing a dy- 
namical instability (see text). 



position x affects the electronic and structural proper- 
ties of the C32 phase of CaA^-^Six. We will consider a 
dense grid of x values in the range [0,2], which includes 
both the physical [0.6, 1.2] and the unphysical range. For 
each x-value considered, the hexagonal unit cell has been 
optimised up to a residual stress of 1 MPa, and band 
structure calculations have been performed in the relaxed 
configuration. For this dense grid of x values a slightly 
lower cutoff of 36 Ryd was used for the plane- wave basis. 
In Fig. we report, as a function of x, the optimised lat- 
tice parameters (upper panel) and g(ep), the calculated 
DOS at the Fermi level (lower panel, black). We start 
our discussion from the experimental x-range. 



Experimental range, 0.6 < x < 1.2 

Taking into account that we have modelled the dis- 
order using the VCA, and that the experiments detect 
several types of impurities, our theoretical in-plane bond 



length appears in reasonable agreement with the experi- 
mental data; the largest discrepancies are observed close 
to x = 0.6, which is the lowest aluminium concentration 
for which the scries has been synthesised. The somewhat 
worse agreement for c is common to the ordered phase 
and may be traced back to a systematic error of the ap- 
proximate exchange-correlation functional. 

In the whole experimental range, a monotonically de- 
creases with x, while c is almost constant. 

The small (less than 3%) variations in the lattice pa- 
rameters correspond to minor modifications in the band 
structure: in the experimental x-range the calculated 
electronic structure of CaA^-zSix (not shown) can be ac- 
curately described by a simple rigid-band model, where 
the Fermi level changes only according to the different 
number of electrons per cell N c (x) = 8 + x, see Fig^ 

An indirect proof of the validity of this model is given 
in the lower panel of Fig. El where we plot, as a function 
of x, the self-consistent DOS at the Fermi level, g(ep), 
of CaA^-ajSia; and the corresponding value for a rigid- 
band filling N c (x) of the (fixed) bands of Ca(AlSi), i.e., 
the bands shown in the middle panel of Fig^ As we can 
see, in the composition range 0.6 < x < 1.2, where the al- 
loy is experimentally found to be stable in this crystalline 
structure, the two curves are almost identical; deviations 
start to appear at the extrema of this range. Then, out- 
side this range, the rigid-band model completely breaks 
down. The dependence of g(ep) on composition thus 
correlates with rigid-band and structural properties. Its 
almost monotonic variation in the experimental range of 
stability 0.6 < x < 1.2 does not, instead, correlate with 
the observed trend in the superconducting transition: the 
experimental T c (a;) has a sharp peak at x = 1 and drops 
off rapidly as the Al/Si ratio is changed in either direc- 
tions^ This suggests, by exclusion, that the experimen- 
tal trend in T c (x) must be associated with changes of 
the phonon frequencies and/or electron-phonon coupling 
strength as a function of composition. 



High aluminium concentration, x < 0.6 

For x < 0.6 the a lattice parameter still increases as x 
decreases, its slope being larger than in the experimental 
range of stability. Much more peculiar is the behaviour of 
the c parameter. Reverting the trend observed in the ex- 
perimental range of stability, c decreases as x decreases, 
until it equals the a value, for x ~ 0.3 — 0.4; then, at 
x = 0.2, it undergoes an abrupt change, such that the 
c/a ratio jumps by ~ 20%. Our lattice dynamical cal- 
culations show that this change is accompanied by the 
occurrence of dynamical instabilities. 

In Fig. 01 we plot, as a function of x, the CaA^-xSLr 
phonon frequencies of the Bi ff branch calculated at T, A 
and L (purely imaginary frequencies are shown as neg- 
ative values). At x < 0.2 there is a large jump in the 
phonon frequencies of the Bi ff branch at A and L, which 
become imaginary (shown as negative), while the zone- 
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center phonons are still all stable. The entire CaAl2 
(x = 0) phonon spectrum, shown in the lower panel of 
Fig- H reveals that the compound is dinamically unstable 
against lattice distorsions with q z — ±tt/c and eigenvec- 
tors corresponding to the B\ g optical and Am acoustical 
branches, both implying an out-of-plane displacement of 
the Al - Si atoms. 



T 




1 1 i I i I i I i I i I i I i I i I i I 

0,2 0,4 0,6 0,8 1 1,2 1,4 1,6 1,8 

x 



FIG. 4: Selected phonon frequencies of the spectrum of 
CaA^-rcSia; as a function of the Silicon content x. The branch 
which corresponds to the out-of-phase displacement of the 
Al/Si sublattice along z (Bi 9 at F, B2 U at A) becomes unsta- 
ble around k z — ±7r/c for x < 0.2 (imaginary frequencies are 
shown as negative). 

The dramatic changes in the lattice parameters and 
phonon frequencies with compositions x < 0.6 can be 
rationalised in terms of the electronic band structure. 
We will start our discussion from the bands calculated 
at x = 0.6, shown in Fig. [3J At this composition both 
the interlayer band and the tt bands cross the Fermi en- 
ergy, while the a bands are still completely full, even 
though, around A, they are now very close to ep- On 
these grounds one may think that, by further decreasing 
x (i.e., by increasing the aluminium content), it would 
be possible to obtain a holes, as in MgB2. This simple 
hint proves, however, incorrect, since, as already men- 
tioned (when discussing the lower panel of Fig. |SJ), the 
rigid band model breaks down for x < 0.6. We are going 
to make this statement more precise. Our first-principles 
calculations have revealed that the a, -k and interlayer 
bands are no longer rigid when the aluminium content is 
increased above 70% (i.e. for x < 0.6). 

The 7r band decreases its binding energy on the k z =0 
plane, progressively losing electrons and getting emptier 
and emptier as x is reduced, while the interlayer and 
the a bands remain partially filled to conserve the total 
number of electrons per unit cell; but then, as x further 
decreases, the crystal structure, and therefore the bands, 
substantially rearrange. In a rigid band picture, based 
on the x = 1 DOS, the interlayer band would become 
completely empty at x=0.5 (see middle panel of Fig^l; 



IDOS (electrons f.u.) 
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DOS (Slates/eV f.u} 

FIG. 5: Energy bands (black), electronic DOS (blue) and 
integrated DOS (solid red) at x =0.6, i.e., the lowest Si con- 
centration for which CaAISi is observed in the C32 structure. 



but, precisely around this composition, the rigid-band 
picture starts to break down, and the behaviour of the 
c/a ratio starts to change in such a way as to keep the 
interlayer band always at 6f- 

When, at x = 0.2, the interlayer band is almost empty, 
our self-consistent calculations predict a further, abrupt 
change in c/a. This sudden relaxation is related to a 
major rearrangement of all bands, whose main result is 
that the interlayer band stays at the Fermi level, while 
the 7T (bonding) band becomes completely empty in the 
k z — plane. 

The complete emptying of the 7r band in the "upper 
half of the Brillouin zone (k z ~ ±7r/c, A — L - H - A 
path) further weakens the bonds in the vertical direction, 
and this explains the appearance of imaginary frequencies 
in the out-of-plane phonon branches of the C32 phase of 
CaA^-^Sia, for x < 0.2. The displacement corresponding 
to these phonons is compatible with the stable structure 
of CaAl2 at ambient condition, the MgCu2 Laves phase, 
in which the Al atoms sit on triangular layers interca- 
lated by mixed Ca-Al planes. In this structure, the main 
bonding force is provided by p — d bonds between Ca and 
Al, while in the C32 structure the in-plane sp 2 bonding 
is important. Interestingly, experiments have showni^ 
that the Laves phase is stable only up to x < 0.2, where 
we find the occurrence of a lattice instability for the C32 
phase. A thorough investigation of the whole phase di- 
agram and transition goes, however, beyond the aim of 
the present study^ 

In the upper panel of Fig. [7] we display the electronic 
bands of CaAl2 corresponding to the optimised C32 struc- 
ture. Because of the large modifications in the lattice 
parameters, these bands bare very little resemblance to 
those of Ca(AlSi), and even less to those of MgB2- What 
seemed to be a good candidate for a compound with prop- 
erties similar to MgB2 turns out to be a very poor im- 
itation, once the lattice parameters are optimised. In 
particular, even though there is a tiny amount of a holes 
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around the A point -the fingerprint and main ingredient 
for superconductivity of MgB 2 -, in CaAl 2 most of the 
DOS at ep comes from the it* bands and from the the 
interlayer band (which in MgB2 is, instead, above ep, 
and thus completely empty). Worst of all, as we already 
saw, CaAl 2 is dynamically unstable in the C32 structure. 

We wish to point out that the differences in the band 
structure and the dynamical instability are both a con- 
sequence of the c/a reduction: if we not only constrain 
the compound to the C32 structure, but also keep the 
c/a ratio at the same value as MgB2 (rather than be- 
ing optimised, as done until now), then we find a holes 
at the Fermi level (see fig. and a dynamically stable 
compound: all the calculated phonon frequencies, not 
shown, are now realrfS 

We finally observe that the abrupt c/a reduction (when 
the crystal is, as in our calculations, artificially forced 
to maintain the C32 structure) amounts to a Fermi-level 
pinning of the bottom of the interlayer band. Since, as 
mentioned, the bottom of the interlayer band has a pre- 
vailing Ca-G? character, this observation in turn suggests 
that the mechanism which makes CaAl 2 so different from 
MgB2 and drives its instability towards the Laves phase 
is, in fact, the pinning at the Fermi level of the calcium 
rf-states. 



High silicon concentration, x > 1.2 

While in the limit of high aluminium concentration 
{x < 0.6) we found several fingerprints and a plausible ex- 
planation for the experimental lattice instability, in the 
opposite limit (high silicon concentration, x > 1.2) the 
only remarkable finding remains the deviation from the 
rigid-band behaviour, which, as x increases, manifests 
itself in the slow descent in energy of the lowest anti- 
bonding a band. This band, which at x = 1 is ~ 1 eV 
above Ep at T, starts filling up at x = 1.6. 

Apart from this deviation (which however is less dra- 
matic than for x<0.6, see lower panel of Fig.HJ), neither 
the optimised structural parameters a and c (upper panel 
of Fig.HJl), nor the phonon frequencies selected in Fig. 0] 
nor, in fact, any frequency of the complete phonon spec- 
trum (not shown), seem to do anything special in the 
composition range 1.2 < x < 2. This suggests that (at 
zero temperature) the C32 phase of CaSi 2 (i.e., x = 2) rep- 
resents a metastable structure, i.e., a local minimum of its 
total energy which is higher than the absolute minimum. 
Such a guess is reinforced by the complexity of the unit 
cells of the actual stable phases of CaSi^, which could 
easily imply sizable energy barriers separating them from 
the metastable C32 phased 

Again, a thorough investigation of the whole phase di- 
agram and transitions of CaSi2 goes beyond the aim of 
the present study. 



Experimental range 



i virtual crystal 
^ c virtual crystal 
a CaAISi 
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a experimental 
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FIG. 6: Upper panel: optimised lattice parameters for the C32 
crystalline structure of CaA^-xSi^. Lower panel: blue cir- 
cles indicate the self-consistent density of states at the Fermi 
level g{cp) of CaA^-^Sia, in the fully relaxed structure, while 
red triangles show g(tF) obtained with a rigid-band filling of 
Ca(AlSi) appropriate to the electron number N e (x) = 8 + 2: 
(see middle panel of Fig. . 



V. SUMMARY AND CONCLUSIONS 

In this work we have reported first-principles calcula- 
tions of the electronic and vibrational properties of the 
ternary alloy CaA^-zSiz in the entire composition range 
< x < 2, discussing how electronic, structural and vi- 
brational properties are related to each other. In this 
context an interlayer band, which is empty in MgB2, 
and whose bottom has Ca-e? character in our compounds, 
plays a major role. 

At x=l the electronic band dispersion and Fermi sur- 
faces, calculated considering both an ordered phase with 
atomic in-plane ordering and a disordered VCA phase, 
confirms the suggestion of Mazin et al&, that the com- 
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FIG. 7: Upper panel: energy bands and DOS for CaAl2 in the 
structurally relaxed C32 structure. Because of the very smalll 
c/a, the 7r bands acquire a large dispersion in the k z direction, 
and are driven well above ef in the k z =0 plane. Notice that 
the resulting Fermi surface (not shown) is still similar to that 
of Ca(AlSi). Lower panel: phonon dispersion and density of 
states corresponding to the same structure. 
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FIG. 8: Band structure of CaAl2 in the C32 structure with 
the same c/a ratio as MgB2 (c/a = 1.142). In this case there 
are a holes at the Fermi level both at A and F and all the 
phonon frequencies (not shown) are real. 



pound crystallises in a completely in-plane ordered phase, 
if synthesised in presence of equal amounts of Al and Si. 
Our calculation of the electron-phonon coupling confirms 
that the soft Bi 5 phonon branch plays a major role in the 
superconductivity of CaAISi, and shows that this mode 
couples the interlayer and 7r* bands. At variance with 
a previous full-potential calculation, we do not find any 
instability in the phonon spectrum of CaAISi. 

At x ^ 1, in the experimentally accessible composi- 
tion range, the CaAla— x Si x electronic bands are accu- 
rately described by a rigid-band model; in this range 
ff(ef) monotonically increases with x, so that the mea- 
sured non-monotonic variations in the superconducting 
properties cannot be associated with changes in the den- 
sity of states at the Fermi level. These findings provide 
a useful reference and a strong motivation for a future, 
complete study of the electron-phonon interaction as a 
function of x. 

The rigid-band model breaks down for x < 0.6, where 
the increasing aluminium content mainly leads to the 
emptying of the n bonding band, which finally drives the 
compound unstable. The occurrence of such a dynami- 
cal instability is related to an abrupt change in the c/a 
ratio and correlates with the persistence of the interlayer 
band at the Fermi level, related to the pinning of the Ca- 
d states. Therefore, although isostructural and isoelec- 
tronic with MgB 2 , the C32 phase of CaAl 2 is unstable 
and doesn't share strong similarities with its electronic 
bands. 

The main result of this work is that the interlayer band, 
empty in MgB 2 , and the "out-of- plane" phonons (optical 
Bi g and acoustical Am at T), irrelevant for supercon- 
ductivity in MgB 2 , play a major role on the stability and 
superconductivity of C32 intermetallic compounds, once 
the structural parameters c and a assume appropriate 
values. This observation goes beyond the CaAISi family 
and represents the starting point for further studies of a 
wider class of hexagonal, graphite-like compounds^ 
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